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Abstract. Surface plasmon-polaritons have recently attracted renewed interest in 
the scientific community for their potential in sub- wavelength optics, light generation 
and non-destructive sensing. Given that they cannot be directly excited by freely 
propagating light due to their intrinsic binding to the metal surface, the light-plasmon 
coupling efficiency becomes of crucial importance for the success of any plasmonic 
device. Here we present a comprehensive study on the modulation (enhancement 
or suppression) of such coupling efficiency by means of one-dimensional surface 
corrugation. Our approach is based on simple wave interference and enables us to 
make quantitative predictions which have been experimentally confirmed at both the 
near infra-red and telecom ranges. 
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1. Introduction 

Surface plasmon-polaritons (SPPs) are electromagnetic modes originated from the 
interaction between light and mobile surface charges, typically the conduction electrons 
in metals [1]. Because of the so-called "excess of momentum" with respect to light 
of the same frequency, SPPs cannot propagate away from a planar surface and are 
thus bound to and guided by it. As a consequence of such binding, SPP modes can be 
laterally confined below the diffraction limit, which has raised the prospect of SPP-based 
photonic circuits [2[ [3], H] . To build up this kind of circuits one would require a variety 
of components in which incident light would be first converted in SPPs, propagating 
and interacting with different devices before being recovered as freely propagating light. 
Hence, a great deal of attention has been recently devoted to the creation of optical 
elements for SPPs El [3, El El [10], as well as to the efficient coupling of freely- 
propagating light into and out of them. This latter issue constitutes the fundamental 
bottleneck that must be overcome in order to fully exploit the potential of SPPs, given 
that established techniques for SPP generation (which make use of prism fH] [12] . 
grating [13] or nanodefect [l^ coupling) require that the system's size be well out 
of the sub-wavelength scale in order to obtain a neat SPP signal. On the other hand, 
p-polarized back-side illumination of sub-wavelength apertures in optically thick metal 
films [l5l[T6llI71[T8lll9l[20l[2ll[22] prevents both damping and signal blinding but it 
does not ensure a unique propagation direction for the generated SPPs. 

In a previous work [23], we proposed a novel back-side slit-illumination method 
based on drilling a periodic array of indentations at one side of the slit. It was 
demonstrated that the SPP beam emerging from the slit to its corrugated side can be 
back-scattered in such a way that it interferes constructively with the one propagating 
in the opposite direction, thus obtaining a localized unidirectional SPP source. Here, 
we provide a comprehensive version of our proposal and discuss in some extent its range 
of validity. Additional experimental measurements will be also presented. 

This paper is organized as follows: in Section [2] we summarize the key concepts of 
our proposal and focus in some quantitative aspects of SPP generation and reflection. 
The validity of our simple wave interference model is discussed in Section [31 Finally, 
experimental results are presented in Section [H prior to general conclusions. 

2. Description of our proposal 

A picture of the proposed structure is shown in figure [H A periodic array of one- 
dimensional (ID) indentations is fabricated at the output metal surface close and parallel 
to the illuminated slit. The starting point for such a design can be found in a previous 
work on ID SPP scattering by means of a modal expansion formalism [Ml [25] . In 
order to cope with SPP launching, we considered a single slit flanked by an array of 
indentations (rectangular grooves) placed in the output surface of a thick metallic film. 
Eventually, the distance between the slit and indentations was taken to be infinity. In 




Figure 1. Scanning electron micrograpii and scliematic diagrams of the proposed 
structure. Parameters {ai,h,w,d,P} defining the geometry of the system are also 
shown. 

this way, the sht merely played the role of a theorist's SPP-launcher, as far as it can be 
shown that the field created by the slit corresponds to SPP illumination into the grooves. 
Besides, we also found a simple geometrical condition for the groove array to behave as 
a perfect Bragg mirror, associated to the low-A edge of the plasmonic band gap for the 
periodic system. Combining these two elements, one can obtain a remarkably simple 
scheme to modulate the SPP coupling-in at a real back-side illumination experiment: 
given an incident wavelength, let us design a groove array for which SPP reflectance 
rises to a maximum and place it at a distance d from the slit (Situation B of figured]). 
Hence, any outgoing SPP generated at the same side of the slit will be mainly back- 
scattered by the grooves and interfere either constructively or destructively with the one 
that is generated at the opposite side. This interference can be tuned by adjusting the 
separation d between the slit and the first groove of the array, defined centre to centre. 
The total phase difference, 0, between the interfering SPPs will then consist of the phase 
change upon reflection plus the additional shift resulting from the two different path 
lengths along the metal: 

(P = ^R + 2Re[kp]d, (1) 

where kp holds for in-plane plasmon wave- vector. According to ([T]), constructive or 
destructive interference should occur for those phase values equal to, respectively, even 
or odd multiples of vr. 

It is clear that, as will be discussed in section [3], several objections may arise against 
this very simplified model, but before we turn to its validity, let us take a closer look at 
the two ingredients on which it is based: the generation of SPPs at a sub-wavelength 
aperture and the phase they acquire as a result of Bragg reflection. 





Figure 2. Calculated |i?e[i/j,]| distribution at the output surface of an Au film 
perforated by a single slit (Situation A of figure [T]) . Incident light is p-polarized and 
impinges normally onto the back side of the metal surface at A = 800 nm . Here, slit 
width oo = 160 nm and film thickness h = 300 nm. 



2.1. SPP generation at a single sub-wavelength slit 

Figure [2] renders a finite-difference time-domain (FDTD) [26] simulation of 
the electromagnetic (EM) field distribution originated from p-polarized back-side 
illumination of a sub-wavelength slit on a thick Au film. As can be seen, the most 
of the field is diffracted away, but a significant fraction appears to be bound to the 
metallic surface at each side of the aperture. However, to what extent is that a confined 
state of surface plasmon? 

Such an assessment requires that analytical expressions for the EM field distribution 
created by the slit be obtained without any a priori assumption about the presence of 
SPPs. For that purpose, we have made use of the above-mentioned modal expansion 
technique. Given that it has been extensively described elsewhere [211 ESj, let us just 
briefly summarize its basic ingredients: the EM fields are expanded in terms of the 
eigenmodes in each spatial region (plane waves at input/output regions and waveguide 
modes inside the indentation) and then the expansion coefficients are obtained by just 
matching appropriately the parallel components of the fields at the two metal-dielectric 
interfaces. The dielectric response of the metal is taken into account by applying surface 
impedance boundary conditions (SIBC) [27] to the tangential components of the EM 
fields at the metallic surface. For a non-magnetic medium, 

Ft{r) = Et{r) - ZsHt{r) x n(r) = 0, (2) 

where Zg = e{\)~^^'^ and n(r) is the unitary vector normal to the surface directed into 
the metal half-space. However, SIBC are not applied at the vertical walls defining the 
slit but for the calculation of propagating constants along the z direction. This choice 
allows us to express the EM fields inside in terms of the waveguide eigenmodes of a 
perfect conductor (PC), which are known analytically. Although the absorption inside 
the cavities is therefore neglected, one can expect this not to be a serious shortcoming 
when considering sizes much greater than the skin depth. The end product of our 
expanding and matching is a linear system of algebraic equations that connect the 
modal amplitudes of the F field at the input and output openings of the slit. Once 
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those self-consistent amplitudes are found, it is straightforward to obtain the EM fields 
at any desired point. 

By imposing the constraint that incident light impinge in "classical mounting" (ie. 
within the xz plane), we just have to concern ourselves with the ^/-component of the 
magnetic field. From a mathematical point of view, Hy at the output side is obtained by 
integrating all across the slit every considered eigenmode 0„ multiplied by a scalar ID 
Green's function and then weighting each contribution with the corresponding amplitude 
E'^ at the output opening: 

rxo+ao/2 

Hy{x, z) = —'^^ E'^ / dx'G{x,x'; z)(j)n{x'), (3) 

^ Jxo-ao/2 

where z stands for the distance from the output surface. This closely resembles the 
Huygens-Fresnel description of wave propagation in terms of a set of punctual emitters, 
but we have to keep in mind that all those "emitters" are self-consistently connected. 

However, the information on the character of the generated field is contained neither 
in the modes nor in their amplitudes, but in the propagator itself: 

where fco = 27r/A. Despite its impressive appearance, G{x,x';z) just computes the 
projection of EM fields at the opening of the slit onto all possible diffracted waves, 
whether they are propagating or evanescent. Bound-to- interface contributions are 
incorporated into the picture as a consequence of finite Zg, which makes the difference 
with respect to perfect conductor approximation: for = 0, (jlj) transforms into an 
integral representation of the Oth-order Hankel function of the first kind (ie. the well- 
known Green's function for two-dimensional Helmholtz operator [28]), otherwise it has 
to be evaluated numerically. Such numerical inspection reveals that G{x, x'; z) tends 
to the PC result for |x — x'l << A irrespective of Zg pH [25]. On the other hand, 
in the regime where z, \x — x'\ ~ 0{X), oscillatory contributions within the kernel of 
(jl]) mutually cancel everywhere but in the region close to the integrand singularities at 
k = ±kp, with kp satisfying 

\lkl-kp^ = -Zsko. (5) 

This is, by the way, the SPP dispersion relation of a flat metal-dielectric interface within 
the SIBC. In that asymptotic limit, the Green's function can be explicitly approximated 

as 

Gas{x,x';z) = -M^e^C^pl-^'l-W^ (6) 
kp 

Therefore, and even in the presence of absorption, SPPs govern the EM coupling along 
the surface at a distance of several wavelengths, whereas "PC-like" behaviour is observed 
at the close vicinity of the slit. It is worth mentioning that this simple fact is completely 
misinterpreted in several recent papers, as pointed out in a previous work [29]. In 
any case, the existence of these two regimes have also been remarked by introducing 
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a "creeping wave contribution" that rapidly vanishes for increasing distances and is 
exphcitly defined as the difference between total and SPP fields along the metal-dielectric 
interface [30]. 

In order to determine the precise range of distances for the EM field at figure [2] 
to be dominated by either "PC-like" or SPP contribution, we have calculated |i?e[ify]| 
at the metal surface for values corresponding to that of Au at 800 and 1500 nm. 
Each calculation was carried out for the exact, asymptotic and PC versions of Green's 
function. As can be seen in figure [3l comparison with the exact result in the near- 
infrared (NIR) shows that the asymptotic limit is already reached for a distance of 
about 2 A from the center of the slit, which is increased up to 6 A when the incident 
wavelength falls within the telecom range. Consequently, it is only for greater distances 
that we can unambiguously establish a one-to-one correspondence between fields at the 
interface and SPPs. 

In figure H] we present the fraction of the output current that is transferred into 
SPPs (/spp) and scattered out of the plane (/out) for the same parameters as in 
figure [3] all across the sub-wavelength regime. Given that SPPs gradually attenuate 
when propagating along the metal, the values for /spp are calculated at x = ±ao/2 
in order to compare with those of /out- As the slit width increases, the out-of-plane 
radiation is clearly favored at the expense of the coupling into SPPs, which can be 
easily found to be proportional to (sin[/cpao/2]//cpao)^ because of the geometry of the 
system |2H |25]. For typical experimental width = 160 nm, no more than 30% of the 
output energy is driven into SPPs at A = 800 nm and such percentage is reduced to 17% 
at A = 1500 nm. These values are in good agreement with those previously reported 
pT] and provide a preliminary estimate of the expected performance for our proposed 
slit+grating structure when operating at perfect constructive interference condition. 

2.2. Phase shift upon Bragg reflection 

As mentioned at the beginning of the section, it has been shown that the refiection 
of SPPs by a periodic array of indentations presents maxima at those frequencies 
corresponding to the low-A edges of plasmonic bandgaps [25] . For narrow sub- 
wavelength indentations, the spectral locations of these edges can be approximated by 
folding the dispersion relation of SPPs for a fiat metal surface into the first Brillouin 
zone [31]. Within the SIBC, such folding results in 

kpP = kQRe[qp]P = mvr, m = 1, 2 . . . (7) 

where P is the period of array and Qp = ^Jl — Z'^. Remarkably, although the refiectance 
maxima depend on the groove geometry (width and depth) and the number of grooves 
[211 [25], their spectral locations do not (see figure [5]). 

Assuming that A and P fulfill ([7]), let us consider the phase shift for a given 
resonant wavelength Ar. Information on such shift is contained in the complex refiection 
coefficient r relating the amplitudes of incident and refiected fields. Although the 
obtention of r is usually regarded as a mere preliminary to that of refiectance (defined 
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Figure 3. Calculated |i?e[ffj^]| as a function of the distance from the center of the 
slit, evaluated at the output surface of an Au film. Solid lines represent the full 
calculation, whereas dashed and dotted ones stand for perfect conductor approximation 
and asymptotic expansion, respectively. The geometrical parameters are: slit width 
flo = 160 nm and film thickness h = 300 nm. Incident light is p-polarized and impinges 
normally onto the back-side of the metal surface. Panels (a) and (b) show results for 
A = 800 nm and A = 1500 nm, respectively. 

as R= we can always establish a straightforward connection between r and phase 
shift (l)R. 

cos 0i? = -Re[r]/|r|; sm(j)R = Im[r]/\r\. (8) 

Once the asymptotic limit is already reached, these auxiliary magnitudes cos(^ij, sin0R 
provide complete information about SPP shift upon reflection, irrespective of the exact 
distance at which fields are evaluated. We have found that is close to vr over a wide 
range of groove depths for a/A < 0.2 at both NIR and telecom ranges, as can be seen in 
figure El Taking this result into account and substituting for kp from ([7]) into (P) yields 

0(Afi) = (2mrf/P + l)7r, (9) 
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Slit width to wavelength ratio (a^/X) 

Figure 4. Fraction of the energy that is transferred into SPPs (lines) and scattered 
out of the plane (symbols) at the output surface of an Au film perforated by a single slit 
that is back-side illuminated with p-polarized light. The values for /spp are calculated 
at distances of ±ao/2 from the center of the slit. 

which reduces the design of our proposed scheme to a suitable choice of the d/ P ratio. 
3. Validity of the simple wave interference model 

Our previous discussion leading to ([9]) implies that slit and grating be considered as 
independent elements. Therefore, it does not take into account the radiation coming 
back from the grooves, while, in principle, EM-fields at all openings have to be self- 
consistently calculated [32]. In order to quantify the "perturbation" of the SPP source 
(ie. the slit), we define a re-illumination parameter ^ that averages the modification of 
the x-component of the electric field inside the slit originated by the adjacent grating: 

i = - dx'\l-E,{x')/E:^{x% (10) 

do J-ao/2 

where qq is the width of the slit, the x-component of the electric field calculated in 
the presence of the array and E^^ the one obtained for the isolated slit. 

In figureHwe present a contour plot of ^ vs groove depth and slit-to-array separation 
for a system with 10 grooves at A = 800 nm. As can be seen, the modification 
of the field pattern within [400, 800] nm is below 15% and ^ rapidly decreases for 
increasing distances, thus supporting our implicit assumption in ((91). With respect 
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Figure 5. Calculated reflectance of SPPs by a finite periodic array consisting of 5 
(solid line), 10 (dashed) and 15 (dotted) grooves carved on Au. Here, a — w — 100 
nm and P = 390 nm. Gray-shaded area marks the region where a plasmonic bandgap 
occurs. SPP fields are evaluated at a; = — 3.5/im (w —S.SXmax), the origin being 
located at the center of the first groove. 



to the dependence on groove depth, it is governed by the reflectance properties of the 
array, ^ rising to its maximum as R does (see inset in figure [7j). Such a maximum 
becomes clearer the more separation approaches to the plasmonic regime {d ~ 3A). 
On the other hand, modulation along the vertical axis results from simple interference 
between counter-propagating SPP waves originated at the slit and its nearest groove. 
Therefore, sequential minima of ^ appear for d = (2m + l)Ap/4, whereas {^max} are 
associated with d = m\p/2, given that Xp = 2t: /kp and m = 0, 1, 2 . . . 

However, the key point of our proposal still relies on SPPs being reflected by a 
groove array, while the EM fields radiated by the slit cannot be considered "purely 
plasmonic" but at a distance of several wavelengths (see figure [3]). In order to 
characterize the efficiency of the sli+array system as an SPP-launcher for any slit-to- 
array separation, we introduce its "efficiency ratio" , Ej^: given that the array be located 
at the left side of the slit (see figure [1]), Er is defined as the quotient between the current 
intensity of right-propagating SPP with and without the grooves. Strictly speaking, Er 
provides the efficiency of the output side of the device. The total efficiency, defined 
as the percentage of incident energy transferred onto the plasmon channel, strongly 
depends on the illuminating setup. Er should vary within the interval [0, 4] showing a 
dependence on the distance between the illuminating slit and the groove array. More 
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Figure 6. Calculated values of cos (j)R{XR), sin0ij(Afl) of a 10-groove array in Au for 
increasing values of groove depth. Panels (a) and (b) display results for a = 100, 150 
and 200 nm (solid, dashed and dotted lines, respectively) evaluated at Xu — 800 nm 
(P = 390 nm, m = 1). Resuhs for a = 200, 300 and 400 nm at = 1500 nm (P = 750 




Groove depth (nm) 

Figure 7. Contour plot of re-illumination parameter ^ in Au as a function of groove 
depth and slit-to-array separation. Here uq = ai = 160 nm, P = 390 nm, iV = 10 and 
A = 800 nm. Vertical dashed line marks the maximum value of the reflectance curve 
for the isolated groove array, which is presented at the lower right inset panel. 
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importantly, Eji > 2 implies that the right-propagating SPP current in the presence of 
grooves is larger than the total SPP current (left- plus right- moving) in the single slit 
case, so some of the power radiated out of plane is redirected onto the SPP channel. 
According to our simple wave interference model, 

Er^ \l + re^'''''^\'^, (11) 

where r is the complex reflection coefficient of the groove array for SPPs. 

To check the validity of iQ and f|TT]) for slit-to-array separations outside the 
asymptotic regime, we have carried out numerical calculations of EM fields by means 
of both modal expansion and FDTD. The system under consideration is intended to 
operate at a wavelength of 800 nm on a gold film [33]. We consider an array of 
10 grooves with a period P = 390 nm. The depth of the grooves is chosen to be 
w = 100 nm, while the width of both grooves and slit is a = 160 nm, which are 
typical experimental parameters. Figure El^a) shows the comparison between ( fTTj) and 
numerical evaluations of Er, as well as the location of interference maxima (vertical 
lines) predicted by ([9]) for m = 1. The agreement between the modal expansion and 
FDTD results is excellent but for distances at which intra-wall coupling between the 
slit and the first groove has to be taken into account {d ^ 2a). As can be seen, the 
locations of maximum Er are accurately predicted by ([9]), which allows us to design SPP- 
launchers without elaborate numerical calculations. Moreover, the simplified model of 
(|TT1) provides a good approximation to Er with the sole input of r. This also implies that 
non-plasmonic contributions to groove illumination play a minor role in the occurrence 
of either constructive or destructive interference, which is clearly described by (fTTI) with 
the except of minor shifts. 

In addition to the efficiency ratio, field patterns in both minimum and maximum 
condition were also calculated using the FDTD method. As shown on figure [H](b), SPPs 
are completely absent from the left side of the slit whereas field intensity at its right 
side is clearly modulated by the slit-to-array separation, which also governs the spatial 
distribution of the field that is radiated into the vacuum. 

4. Experimental results 

4.I. Near-infrared measurements 

For our proposal to be tested out at the NIR regime, several slit+array samples were 
fabricated on gold films with a Focused Ion Beam (FIB). As described in [23], each 
sample consists of a single long [L = 30/im) slit of width a^ = 160 nm perforated at a 
300-nm-thick film that is fianked by a periodic array of grooves (P = 390 nm, a = 160 
nm, w = 100 nm). Such array is placed at a given distance d and only extends over L/2 
(see figure [I]). This kind of samples enables us to measure Er, as the upper part can be 
used as an in-chip reference of the "isolated slit" . 

A set of samples with d = {195,292,390,486,585} nm was imaged at 800 nm by 
a Photon Scanning Tunneling Microscope (PSTM) making use of an incident focused 
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Figure 8. Numerical results for the SPP launcher at wavelength A — 800 nm. (a) 
Dependence of the efficiency ratio Eji on the slit-to-array distance. The geometrical 
parameters defining the system are: slit and groove widths a = 160 nm, groove depth 
w = 100 nm and array period P = 390nm. The figure renders the curves obtained by 
means of FDTD (solid), modal expansion (dashed) and equation (jlip (short-dotted). 
Vertical lines mark the positions of En maxima according to ([9]). (b) Calculated 
|i?e[7?j,]| distributions over xz plane for two different distances corresponding to 
minimum and maximum values of at A = 800 nm. 
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beam illumination. For each sample, a pair of images was recorded by scanning at a 
constant distance of about 60 to 80 nm from the surface (see figure [9](a)). The first 
image of the pair, corresponding to a SPP generated from a single slit, is obtained by 
focusing the laser beam on the upper part of the slit. For the second image, the laser 
beam is moved to the lower part in order to collect the data for the slit+array structure. 
An average longitudinal cross-cut of each image is obtained by using 20 longitudinal 
cross-cuts, corresponding to different coordinates along the slit axis. Then, the relative 
position of the two average cross-cuts is adjusted so that the saturated areas (i.e. the 
signal taken right on top of the slit) are super-imposed. Finally, the experimental 
efficiency ratio, Efj, is extracted by averaging the ratio between the two curves along 
the longitudinal cross-cut. Figure [9)d renders experimental values (circles) of for the 
five different samples fabricated, as well as the ones obtained from FDTD simulations 
(solid line). The concordance between measurements and theoretical predictions is 
quite remarkable, especially when taking into account that each experimental point 
corresponds to an average over a different set of samples. We find that this agreement 
(previously reported in [23]) provides a clear support to our proposal for a localized 
unidirectional SPP source. 

Another way of looking at the role of surface corrugation is to consider its infiuence 
on the fraction of the output energy that is radiated into vacuum. Given that some of 
the radiated power is redirected onto the SPP channel for the condition of maximum 
Er (see field pattern at figure [8]), we may wonder whether or not the radiated field is 
also modulated by the sht-to-array separation. For that purpose, a new magnitude Eout 
can be defined as the ratio between the radiated energy with and without the grooves. 
According to our numerical simulations, such "out-of-plane efficiency" presents a similar 
(but opposite) dependence on d to that of Eji. In order to obtain experimental values 
for Eout, a new type of sample was designed (see figure [TOTa)). Now, the illuminating 
slit is fianked by two groove arrays with the same periodicity P = 390 nm, each one 
extending over L/3. No corrugation is present at the middle part of the system, for it 
to be used as the "single slit" reference. The array on the top is located at a distance 
di = 607 nm for which the coupling to SPPs rises to a maximum at A = 800 nm, whereas 
a minimum appears for the distance d2 = 404 nm of the bottom one. Consequently, the 
far-field radiation pattern of the composed structure is expected to present a c?i — ^2 
ascending staircase profile. 

In figure fTOlfb). (c) we present PSTM images recorded at 800 and 3000 nm from the 
surface of the sample. As can be seen, the intensity distribution along the illuminating 
slit increases from the upper to the middle third, as well as from the middle to the lower. 
Although this behavior is in qualitative agreement with our predictions, a rigorous 
determination of E^ut would have required extensive measurements similar to those 
of Ef{. Unfortunately, such procedure became impossible because of accidental fatal 
damage in the sample. However, we have managed to obtain a rough estimate of Eout 
from available PSTM images: dashed rectangles in figure fTOT b). (c) mark the areas over 
a 4-line average (~ 625nm) longitudinal cross-cut of each image is obtained by means 
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Figure 9. Experimental measurement of at A = 800 nm for the same geometrical 
parameters as in figure [H (a) PSTM micrographs recorded for a sample with d = 585 
nm at both "single slit" (top) and slit+array configurations (bottom). The right panel 
shows the two cross-cuts from which Ef{ is obtained. Vertical lines define the interval 
along the ratio is averaged, (b) Experimental (circles) and numerical (solid line) values 
of En as a function of slit-to-array distance. The error bars represent the standard 
deviation over a set of different structures with the same nominal parameters. 



of WSxM software [34]. The resulting intensity profiles at figure [TOTd). (e) show a clear 
succession of steps, which we decide to characterize by the arithmetic mean along the 5- 
micron central segment of each plateau. Numerical estimates of Eout are then calculated 
as the ratio between di, d2 and single slit values (see table at the bottom left of figure fTOl). 
The coincidence of those estimates with the calculated Eout curve in figure [TOT f) is 
amazingly good, which encourages us to carry out conclusive measurements in the near 
future. With respect to figure [TOlff). we finally have to remark that the radiative-to-SPP 
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conversion seems to be more efficient than its opposite, as far as Eout < 2 for any d. 
4-2. Telecom measurements 

Similar samples to those used in the NIR measurements were designed to operate at the 
telecom range by up-scaling the period of the array and its separation from the slit (see 
figure fTTlfa)). However, in this wavelength regime, we found a kind of instability in the 
illumination setup that resulted in a noticeable variation of SPP intensity during the 
near-field scan process, which takes about 45 minutes per image. As a consequence of 
those intensity jumps, the technique used to evaluate the "efficiency ratio" in the NIR 
became unsuitable. Instead, we found Eji as the SPP signal ratio taken from each pair 
of near-field images (with and without side grooves) at the same distance from the slit, 
where its non-plasmonic field contribution can be disregarded, whereas the SPP signal 
is still substantial for the quantification (^ 50/im). To decrease the uncertainty of thus 
obtained efficiency, a series of scans were performed for every structure and wavelength 
measurements, conducting independent adjustments, with the subsequent averaging of 
the En values obtained. Hence, the error of Er represents a statistically estimated 
deviation. 

A typical pair of near-field optical images is presented in figure [TTlfb). (c). For 
telecom wavelengths, the SPP propagation length is increased up to ~ 200/im. Panel 
(c) features a strong SPP beam propagating away from the slit in the direction opposite 
to the array and thereby demonstrating unidirectional SPP excitation. Averaged results 
and estimated errors for Eji (previously reported in [23]) are rendered in figure [TTT d). 
Notice that the validity of our proposal is now tested in a different way: for a given 
slit-to-array separation, Eji is measured within the wavelength range 1500-1620 nm, so 
that the phase difference described by ([1]) is changed with the increasing wavelength, 
providing the conditions for constructive or destructive interference. Obviously, this 
spectral dependence of the efficiency is different for different slit-to-array separation, and 
we support that experimentally. For the case of the sample with d = P + P/2 = 1125 
nm, Er decreases as the wavelength increases (with the only exception of a sharp peak 
at 1520 nm), evolving from a favorable regime {Eji ^ 2) to one in which coupling into 
SPPs is clearly diminished by the array ( Efi < 1). Conversely, Eji ^ 2 for the sample 
with d = 3P/4 = 562 nm all over the range. As can be seen, the comparison between 
experiments and modal expansion calculation is rather satisfactory. 

Finally, we have to mention that the proposed approach for the excitation of 
localized unidirectional SPP beams can also be combined with the appropriate design 
modifications to create functional components for SPP focusing to a spot or tuning the 
SPP beam divergence. If E^i > 2 is expected for a given slit-|-array set, its circular 
bending may produce a converging gaussian beam whose waist length and radius can be 
adjusted by means of the curvature. Several curved SPP focusers has been previously 
achieved [35], IMl ETj [38l |33|, but we find the mirror-blocked back-propagation to be a 
plus. Although the rigorous modelling of SPP coupling at curved structures is rather 





Figure 10. Experimental estimate of Eout at A = 800 mii. (a) Scanning electron 
micrograph of the sample. The geometrical parameters are: slit length L — 30/xm, slit 
width flo — 104 nm, groove width a = 75 nm, groove depth w = 100 nm and array 
period P — 390 nm. (b), (c) PSTM micrograph recorded at distances of 800, 3000 
nm from the metal surface, (d), (e) Average longitudinal cross-cuts along dashed 
rectangles in images (b) and (c). Vertical dashed lines mark the position of the slit, 
(f ) Calculated i?ojit as a function of slit-to-array distance. Lower left table: summary 
of the experimental estimates of Eout- 
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Figure 11. Spectral dependence of Eji at the telecom range, (a) Scanning electron 
micrograph of the sample. The geometrical parameters are: slit length L = 50/im, slit 
width ao = 400 nm, groove width a = 200 nm, groove depth w = 100 nm and array 
period P = 750 nm. (b) Near field image recorded with the laser beam focused at the 
"isolated slit" position of a sample with d = 3P/2 = 562 nm. (Size = 70 x 2Qfim^, 
A = 1520 nm). (c) Same for slit+array focusing, (d) Spectral dependence of Efi for 
slit-to-array distances oi d = 3P/2 = 1125 nm (experiment: squares; theory: solid 
line) and d — 3P/4 = 562 nm (experiment: circles; theory : dashed line). 





Figure 12. (a) Scanning electron micrograph of the curved structure, characterized 
by sht and groove widths of 400 and 200 nm, respectively, groove periodicity P = 750 
nm, groove depth w = 100 nm and sht-groove distance d — 1125 nm. Film thickness 
h = 280 nm, curvature radius R = 30/im and slit chord length L = AOfim. (b), (c) 
Same for R — Adfim and R = 60/im. (d), (e), (f) Far field images of SPPs excited 
on the structures (a), (b) and (c), respectively, recorded with a charge-coupled device 
camera. 

complicated and falls out of the scope of the present work, we expect ([9]) to still provide 
a good estimation for the proper design of the structure, at least as a starting point. 
On that assumption, we have fabricated several samples consisting of an arc-of-a-circle 
slit flanked by the corresponding array of parallel bent grooves (see figure [T^ (a)-(c)). 
Geometrical parameters aQ,a,w,P are the same as in figure [TT| whereas slit-to- array 
distance is set to d = 3P/2 = 1125 nm. 

As shown in figure fT2l(d)-(f). the effect of SPP launching and focusing can be 
appreciated already at the stage of far-field adjustment due to weak out-of-plane SPP 
scattering by surface roughness. Near-field images of SPP excitation on those structures 
recorded at free-space wavelength of 1520 nm are presented in figure [131 These images 
clearly demonstrate the property of a curved slit to excite convergent SPP beam, with 
the effect being sufficiently enhanced due to the side grooves (cf. [36], E?]). With the 
smallest radius of curvature (30 /im), focusing to a confined spot having size 3x3 
/im^ is observed (see the cross cuts in the lower left panel of figure fT3i) . The SPP beams 
excited on the less curved structures feature an extended waist (figure [T3l(b). (c)), which 
scales (at least visually) according to expectations, providing a wider, and hence less 
divergent, SPP beam. That might be useful for particular applications, e.g. in sensing 
of elongated biological samples or in coupling to low-numerical-aperture waveguides. 
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Figure 13. (a), (b), (c) Near-field images (size 64 x 32/im^) of SPPs excited on tlie 
structures in figure [T2l at A ~ 1520 nm. Lower left panel depicts cross cuts obtained 
from (a) by dissecting the SPP focal spot along longitudinal and transversal directions. 



5. Conclusions 

In conclusion, we have studied the SPP coupling-in at sub-wavelength apertures 
with back-side illumination, presenting a novel proposal for the modulation of such 
coupling-in by means of a finite array of grooves. Our approach is based on a simple 
wave interference model that, irrespective of the simplified description of some of the 
physics involved, has been found in good agreement with both sophisticated computer 
simulations and experimental measurements at NIR and telecom ranges. We find this 
to constitute a stimulating challenge for further developments on a wide range of SPP 
devices. 
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